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ABSTRACT 


The  AFOSR  piogram  at  UCSD  contains  mainly  two  areas  of  optical  information  processing 
research: 

(i)  Optical-analog/electronic  hybrid  computing 

In  hybrid  optical- analog/electronic  compuring^'Wc  hay^  performed  several  new  studies  for 
optical  image  processing,  optical  pattern  recognition  and  optical  solution  of  partial  differential 
equations.  For  example.in^mage  processing  we^esigned  and  fabricated  space  variant  filters  *o 
optically  implement  numerous  space-variant  transformations  (e.g..  Hough  transform  for 
detection  of  high-order  parametric  curves,  coordinate  transforms  to  detect  rotation  and  scale 
invariant  features  of  images,  etc.).  For  spike  removal  from  noisy  images,  wc-hay^  developed  a 
new  parallel  algorithm  suitable  for  optical-analog/electronic  hybrid  implementation.  This 
parallel  algorithm  is  space-variant  and  performs  better  than  the  space-invariant  low-pass  filter 
and  the  time-consuming  median  filter.  We  also  ^mdiqd^  a  single  beam  interferometer,  with  a 
special  photorefractive  BaTiO^  crystal,  which  effectively  subtracts  an  exponentially  weighted 
history  of  the  input  from  the  current  value,  thus  functioning  as  a  time-differentiating  (or  novelty) 
filter.  The  single  beam  interferometer  uses  signal  beam  depletion  from  photorefractive 
amplification  of  noise  (fanout)  in  a  specially  designed  and  cut  BaTiO^  photorefractive  crystal 
(with  input  faces  cut  parallel  to  (100),  (01 1)  and  (01 1)  crystallographic  planes). 

For  optical  pattern  recognition  we  have  completed  the  comparison  study  of  statistical 
pattern  recognition  algorithms  by  investigating  the  performances  of  their  discriminant  filters 
when  implemented  by  phase  only  (PO)  filters.  We  found  that  for  each  of  the  pattern  recognition 
algorithms  studied  (e.g.,  least-square  linear  mapping  technique  (LSLMT)  simplified  LSLMT 
(SLSLMTa  Hotelling  Trace  Criterion  (HTC),  Foley-Sammon  (FS)  transform,  and  Fukunaga- 
Koontz  (FK)  transform)  the  peak  of  the  correlation  intensity  increases  and  the  size  of  the 
correlation  spot  decreases  when  a  PO  filter  is  used.  The  signal-to- noise  ratio  was  improved  for 
the  LSLMT,  SLSLMT  and  the  FK  filters,  but  deteriorated  for  the  HTC  and  the  FS  filters. 
Continuing  our  search  for  new,  more  robust  optical  pattern  recognition  algorithms,  we  have 
formulated  a  nonlinear  transformation  based  on  minimization  of  a  special  energy  function,  (i.e., 
self-organization  neural  network).  In  contrast  to  those  supervised  learning  algorithms  based  on 
data  of  known  statistics,  this  unsupervised  learning  algorithm  does  not  require  a  priori 
knowledge  of  the  statistics  of  the  input  data.  To  implement  the  nonlinear  transformation  for 
optical  pattern  recognition  we  are  constructing  a  hybrid  optical/electronic  system  based  on 
random-phase  coded  matrix-tensor  multiplier,  which  can  provide  in  parallel  and  in  real  time  a 
large  number  (e.g.  -10^  - 10^)  of  adaptive  optical  connections  and  therefore  allow  operation  on 
large  size  images. 
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For  the  optical  solution  of  partial  differential  equations  by  a  confocal  Fabry-Perot  system  with 
optical  feedback  and  photorefractive  gain,  we  have  shown  that  the  analog  accuracy  of  the  optical 
solutions  to  panial  differential  equations  can  be  improved  by  implementing  an  iterative  error 
correcting  algorithm. 


(ii)  Digital  optical  computing  and  nonlinear  optical  devices: 

In  digital  optical  computing  research,  we^tave^  been,  developm^ nonlinear  optical  devices 


and  studying,  their  uses  in  parallel  architectures  for  implementation  of  parallel  algorithms.  The 
nonlinear  optical  devices  are  based  on  integrating  silicon  with  PLZT.  tlsiTl^LPCVD  techniques 


to  deposit  polycrystalline  Si  onto  the  electro-optic  PLZT  substrate  and  then  recrystallizing  the 
polysilicon  using  an  Ar*  laser^  to  produce  single  crystal  silicon  grains,  we  fabricated  and 
tested  a  12  X  12  electrically  addressed  spatial  light  modulator  (SLM)  array  and  a  16  x  16 


optically  addressed  SLM  array.  •)  .5 . 


To  study  optical  parallel  computing  architectures^  w^investigated  different  interconnection 
topologies  and  have  generate^  holographic  optical  elements  to  implement  fully  interconnected 
and  hypercube  interconnected  topologies.  , 
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1.  Introduction 


The  AFOSR  program  at  UCSD  consists  of  two  areas  of  optical  information  processing:  (i) 
optical-analog,/electronic  hybrid  computing  for  optical  image  processing,  optical  pattern 
recognition  and  solution  of  partial  differential  equations,  and  (ii)  nonlinear  optical  devices  and 
digital  optical  computing.  The  optical-analog/electronic  digital  hybrid  system  consists  of  an 
analog  optical  processor  with  input/output  interfaces  to  a  microcomputer.  The  analog  optical 
processor  is  employed  to  perform  time  consuming  computations,  while  the  logical  decisions  and 
controls  are  provided  by  an  electronic  microcomputer.  For  image  processing  we  have 
implemented,  with  hybrid  systems,  numerous  space-variant  transformations  (e.g..  Hough 
transforms  for  detection  of  high-order  parametric  curves,  coordinate  transforms  for  rotation  and 
scale  invariant  feature  extraction,  etc.).  For  spike  removal  from  noisy  images  we  have  deveiOt>cJ 
a  new  parallel  algorithm  for  optical  implementation,  whose  performance  is  superior  to  that  of 
low-pass  or  median  filters.  For  motion  detection,  image  clutter  removal  and  novelty  filtering  we 
have  constructed  and  experimentally  evaluated  a  time-differentiating  (or  novelty)  filter,  which 
consists  of  a  single  beam  interferometer  based  on  a  specially  designed  BaTiOj,  photorefractive 
crystal  (with  input  faces  cut  pirallel  to  (100),  (Oil)  and  (Oil)  crystallographic  planes).  For 
optical  statistical  pattern  recognition  we  studied  the  implementation  of  discriminant  functions  of 
our  phase-coded  processor  using  phase  only  filters.  We  also  formulated  a  new  nonlinear 
transform  for  unsupervised  optical  pattern  recognition,  which  is  presently  being  implemented  by 
a  hybrid  optical/electronic  system  based  on  a  random-phase  coded  matrix-tensor  multiplier. 
Finally,  for  the  optical  solution  of  partial  differential  equations  with  the  hybrid  system  based  on 
a  confocal  Fabry-Perot  processor,  we  have  shown  that  the  analog  accuracy  of  the  optical 
solutions  can  be  improved  by  implementing  an  iterative  error  correcting  algorithm. 

In  the  digital  optical  computing  research,  we  have  been  developing  nonlinear  optical 
devices  and  studying  their  uses  in  parallel  .architectures  for  implementations  parallel  algorithms. 
Nonlinear  optical  devices  based  on  integrating  silicon  with  PLZT  were  fabricated  and  tested 
(e.g.,  a  12  X  12  electrically  addressed  spatial  light  modulator  and  a  16  x  16  optically  addressed 
spatial  light  modulator).  In  cur  attempt  to  construct  digital  optical  processing  systems  we  have 
fabricated  optical  holograms  for  implementation  of  architectures  of  fully  and  hypercube 
interconnection  topology.  Now,  we  are  investigating  the  replacement  of  optical  holograms  by 
computer  generated  holograms  for  programmable  optical  interconnects  and  volume  holograms 
in  photorefractive  crystals  for  reprogrammable  optical  interconnects. 
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2.  Hybrid  (Optical-Analog/Electronic  Micro-)  Computer 

2.1.  Optical  Image  Processing 

2.1.1.  Optical  space- variant  filters  for  image  transformations 

In  the  study  of  digital  image  processing  algorithms,  which  would  benefit  from  optical 
implementation,  we  have  constructed  a  computerized  optical  system  for  design  and  generation  of 
space- variant  holographic  filters  (see  Fig.  1)  [1].  Employing  this  system,  space-variant  filters 
containing  256  x  256  sub-holograms  have  been  generated  in  649F  and  dichromatic  gelatin 
plates  for  optical  Hough  transforms  of  high-order  parametric  curves  and  for  coordinate 
transforms.  These  filters  have  been  used  in  the  real  time  space-variant  optical/electronic 
processor  of  Fig.  2  to  detect  the  three  parameters  of  a  circle  and  the  four  parameters  of  an 
ellipse;  the  experimental  performance  is  summarized  in  Fig.  3  and  Fig.  4  respectively. 
Furthermore,  to  perform  rotation  and  scale  invariant  image  transformation,  we  have  generated  a 
In  r-6  coordinate  transform  filter  for  images  of  size  256  x  256  pixels.  The  experimental  results 
on  performance  of  this  filter  are  shown  in  Fig.  5  and  Fig.  6.  Details  can  be  found  in  Ref.  1. 

2.1.2.  Parallel  optical  algorithm  for  image  spike  noise  removal 

We  have  developed  a  new  parallel  algorithm  for  image  spike  noise  removal  suitable  for 
optical  processing  implementation.  The  new  algorithm  is  space-variant  and  performs  better  than 
the  simple  space-invariant  low-pass  filter  used  to  reduce  noise  in  low-level  image  processing, 
(see  Figure  7(b)),  and  the  time-consuming  median  filter  (see  Figure  7(c)).  The  space-variant 
algorithm  for  spike  noise  removal  can  be  applied  to  a  large  size  image  and  is  based  on 
comparing  the  intensity  of  a  pixel  /(/,;')  with  the  average  intensity  of  its  nearest  neighbors 
(/(/,;')).  If  \I  (i,y)-/(/,y)  (  is  greater  than  the  local  variance  of  the  nearest  neighbors,  then  we 
replace  I  {ij)  by  /(/,;),  otherwise  the  original  value  of  I  (ij)  is  preserved.  Computer  simulation 
results  (see  Figure  7(d))  show  that  the  proposed  algorithm  removes  spike  noise  without  blurring 
the  image.  Table  I  summarizes  quantitatively  the  performance  of  the  proposed  algorithm  and 
compares  it  with  the  conventionally  used  spike  noise  removal  algorithms  (i.e.,  median  and  low 
pass  filter). 

2.1.3.  Optical  time-differentiating  filter  for  image  motion  analysis 

We  have  demonstrated  a  single  beam  interferometer  which  effectively  subtracts  an 
exponentially  weighted  history  of  the  input  from  the  current  value,  thus  functioning  as  a  novelty 
filter  [2].  TTie  filter  uses  signal  depletion  from  photorefractive  amplification  of  noise  (fanout). 
Input  which  is  stationary  for  longer  than  the  photorefractive  response  time  is  strongly  depleted, 
but  any  portions  of  the  input  which  change  their  phase,  amplitude,  wavelength  or  polarization 
are  immediately  transmitted.  The  gain  material  used  was  a  crystal  of  photorefractive  barium 
titanate  with  its  surfaces  on  the  (1,0,0),  (0,1,1)  and  (0,1,1)  crystallogiaphic  planes,  so  that  the 
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optic  axis  lies  at  45  degrees  to  the  input  face,  allowing  access  to  the  peak  gain  with  normal 
signal  incidence  [3].  This  crystal  provided  higher  gain  and  faster  response  than  the  normally  cut 
material.  To  examine  the  time-differentiating  filter’s  real  time  operation,  we  used  a  hybrid 
optical/electronic  system  with  input/output  interfaces:  a  liquid  crystal  light  valve  is  used  to 
convert  a  video  image  into  a  phase  and/or  amplitude  modulated  input  signal  and  a  CCD  camera 
is  used  as  an  optical  to  electronic  converter.  The  experimental  results  are  shown  in  Fig.  8.  A 
contrast  ratio  of  100  was  obser\ed  between  moving  and  stationary  images,  with  a  resolution 
limited  only  by  that  of  the  light  valve.  Potential  applications  of  this  filter  include  image  clutter 
removal,  motion  detection  and  analysis,  tracking,  edge  enhancement,  and  image  time 
differentiation. 

2.2.  Optical  Pattern  Recognition 

2.2.1.  Comparative  Study  of  Phase  Only  Filter  for  Statistical  Pattern  Recognition 

We  have  recently  completed  a  study  of  phase-only-filtering  applied  to  discriminant 
functions  provided  by  various  statistical  pattern  recognition  algorithms.  Phase-only  (PO) 
filtering  employs  only  the  phase  of  the  filtering  function  while  the  amplitude  is  set  equal  to  1  at 
all  points  on  the  filter.  We  performed  computer  simulations  in  order  to  evaluate  the  effect  of  PO 
filtering  as  compared  to  amplitude  and  phase  (AP)  filtering.  The  performances  of  PO  and  AP 
filtering  were  evaluated  using  the  following  characteristics:  peak  correlation  intensity, 
correlation  spot  size,  signal-to-noise  ratio  and  optical  efficiency.  We  were  interested  in 
measuring  these  characteristics  for  the  discriminant  functions  provided  by  various  statistical 
pattern  recognition  algorithms  that  have  been  under  investigation  at  UCSD.  The  five  different 
algorithms  we  evaluated  in  conjunction  with  the  PO  filtering  study  were:  the  Least  Squares 
Linear  Mapping  Technique  (LSLMT)  [4],  the  Simplified  Least  Squares  Linear  Mapping 
Technique  (SLSLMT)  [5],  the  Hotelling  Trace  Criterion  (HTC)  [6],  the  Foley-Sammon 
Transform  (FS)  [7]  and  the  Fukunaga  Koontz  transform  (FK)  [8].  The  training  and  test  data  we 
used  consisted  of  grey  tone  images  of  64  x  64  pixels  with  one  byte  of  grey  tone  data  per  pixel. 
Our  training  set  was  composed  of  5  such  images  of  tanks  and  5  of  trucks.  Our  test  data  set  was 
composed  of  4  tank  and  4  truck  images. 

For  each  algorithm  we  created  two  filters  (i.e.,  PO  and  AP)  using  the  two  classes  of  training 
images.  The  test  images  were  then  used  to  compute  the  various  evaluation  characteristics.  The 
simulation  results  are  summarized  in  Tables  2-6  for  evaluation  of  LSLMT,  SLSLMT,  HTC,  FS 
and  FK  transforms.  The  results  obtained  indicate  that  for  all  algorithms  the  peak  correlation 
intensity  increases  and  the  correlation  spot  size  decreases  when  a  PO  filter  is  used.  The  signal- 
to-noise  ratio  improved  for  the  LSLMT,  SLSLMT  and  the  FK  filters  but  deteriorated  for  the 
HTC  and  the  FS  filters.  The  optical  efficiency  did  not  show  a  significant  change  in  either 
direction. 
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2.2.2.  Nonlinear  Transformations  for  Pattern  Recognition 

The  feature  extraction  process  of  a  pattern  recognition  system  car.  be  considered  as  a 
function  /  that  maps  a  set  of  input  patterns  x  onto  a  set  of  features  /  (x).  The  criteria  usually 
used  for  designing  the  feature  extraction  function  is  to  minimize  an  error  energy  function  E, 
which  is  defined  as  11/  (jc)  — r(A:)ll^,  where  r(x)  is  the  desired  output  feature. 

Conventional  pattern  recognition  systems  [4-9]  use  pre -designed  and  fixed  linear  filters  /, 
which  are  implemented  based  on  a  priori  knowledge  of  the  input  patterns.  If  an  adaptive  pattern 
recognition  system  is  desired,  the  function  /  is  implemented  with  a  recursive  algorithm  [10-12], 
where  an  on-line  learning  rule  is  enforced  to  adjust  /  successively.  Usually  the  function  r(x)  is 
determined  by  the  statistics  of  the  input  pattern.  Therefore,  for  data  of  known  probability 
density  function,  r(x)  can  be  considered  given  and  the  adaptive  process  can  be  seen  as 
supervised  learning.  However,  those  systems  will  not  work  well  in  a  changing  environment, 
where  the  statistics  of  the  input  patterns  are  usually  not  predictable,  or  do  not  obey  the  simplified 
gaussian  assumption.  Furthermore,  for  some  complicated  patterns,  due  to  nonlinear 
characteristics,  patterns  may  not  be  well  classified  by  the  most  frequently  used  linear 
disenminant  function  in  the  feature  domain.  Therefore,  nonlinear  transformations  or  nonlinear 
feature  extraction  processes  should  be  employed  to  extract  more  effective  features  in  order  to 
gain  higher  accuracy  and  reliability  in  pattern  classification. 

We  have  already  shown  that  an  adaptive  pattern  recognition  system  may  be  put  in  terms  of 
neural  network  systems  as  an  adaptive  hetero-associative  memory  [12].  The  learning  process  of 
a  hetero-associative  memory  is  equivalent  to  the  supervised  on-line  training  of  an  adaptive 
pattern  recegr.itior  syc'em,  while  rec..IIing  infermatien  from  a  hctcro  associative  mem.ory  is 
equivalent  to  the  feature  extraction  process  of  an  adaptive  pattern  recognition  system.  Both  the 
adaption  rule  of  a  pattern  recognition  system  and  the  learning  rule  of  a  hetero-associative 
memory  are  methods  that  minimize  the  errors  between  the  actual  outputs  and  the  desired 
outputs. 

To  overcome  the  problems  of  conventional  pattern  recognition  systems  as  stated  above,  we 
have  been  developing  a  neural  network  which  utilizes  the  concept  of  self-organization  [1 1,13,14] 
to  optimize  the  solution  without  having  input  statistics.  The  proposed  neural  network  is  based 
on  the  minimization  of  a  special  energy  function,  £' h ll/2(/i(x))-x|l^,  where  both/i  and/2 
are,  in  general,  nonlinear.  In  fact,  the  energy  function  is  a  mean-square  difference  (error) 
between  the  true  pattern  x  and  the  estimate  pattern  x  =  Without  knowing  any  input 

statistics,  the  network  will  self-organize  its  interconnections  by  adaptively  modifying  their 
strength  in  a  direction  that  minimizes  the  energy  function  E'.  This  energy  minimization  process 
will  result  in  optimum  transformations  f\  and  /2,  so  that  the  resultant  self-organizing  feature 
vector  f\{x)  will  optimally  describe  the  input  data  statistics.  This  is  different  from  the  self¬ 
organization  algorithm  developed  from  the  gradient  method  [11,13]  and  the  algorithm  based  on 
cooperative-competitive  learning  [14]. 


The  energy  minimization  process  for  self-organization  will  require  the  development  of 
techniques  capable  of  finding  the  desired  optimum  nonlinear  transformations  f\  and  /2.  In 
supervised  learning  one  of  the  widely  used  techniques  is  the  delta  rule  [9].  We  are  studying  the 
extension  of  the  delta  rule  as  applied  to  the  minimization  of  the  energy  function  £’  for 
unsupervised  learning.  Another  technique  used  in  unsupervised  learning  is  the  simulated 
annealing  algorithm  which  is,  in  fact,  a  stochastic  approach.  The  simulated  annealing  algorithm 
is  well  suited  for  the  minimization  of  the  proposed  energy  function,  and  therefore  is  also  under 
investigation. 

To  implement  the  discussed  algorithm  optically  we  are  presently  constructing  an  optical 
neural  network  architecture  based  on  a  real  time  random  phase  coded  matrix-tensor  multiplier 
interfaced  to  a  microcomputer.  This  optical  memory  (see  Fig.  9)  will  be  capable  of  storing  a  very 
large  number  of  optical  connections  (-10^-10^)  using  a  photoref* active  crystal.  All  these 
connections  can  be  accessed  in  parallel  and  in  real  time  to  provide  the  output  matrix  for 
electronic  thresholding  and  calculation  of  a  correction  term  for  the  tensor  according  to  the 
learning  rule.  The  correction  term  will  be  used  to  update  the  optical  interconnection  stored  in 
the  photoreffactive  memory  material. 

2.3.  Solution  of  Partial  Differential  Equations 

The  Bimodal  Error  Reduction  (BER)  algorithm  is  used  for  improving  the  accuracy  of  the 
solutions  obtained  from  a  moderate-accuracy  processor  [15,16].  Our  investigation  involves  the 
experimental  application  of  this  algorithm  to  satisfy  the  desired  accuracy  requirements  on  an 
electronic/optical  hybrid  processor  being  constructed.  The  optical  portion  of  this  hybrid 
processor  performs  the  analog  solution  of  Poisson’s  equation  for  a  two-dimensional  input  image. 
The  electronic  ponion  updates  the  input  image  and  performs  support  functions.  We  have 
theoretically  analyzed  the  error  correction  algorithm  applied  to  Poisson’s  operator,  which  led  to 
an  expression  for  the  condition  of  convergence.  Also,  a  computer  simulation  of  the  hybrid 
(elcctronic/optical)  implementation  of  the  algorithm  wa^  performed  to  confirm  the  theoretical 
analysis.  The  theoretical  development  considered  additive  error  on  the  operator  in  the  analog 
processor  and  on  the  input  image  presented  to  the  analog  processor.  Both  error  functions  were 
assumed  to  be  spatially  and  temporally  constant  in  this  development.  Conditions  for 
convergence  of  the  algorithm  under  these  assumptions  were  determined  and  found  to  include  the 
expected  operating  conditions  of  the  hybrid  (analog/digital)  processor  being  constructed. 

A  simulation  of  the  BER  algorithm  as  implemented  on  the  hybrid  processor  confirmed  the 
theoretical  results  (see  Table  7).  Refinement  of  the  solution  to  8-bit  accuracy  is  obtained  in  four 
iterations.  Next  we  enhanced  the  computer  simulation  to  more  realistically  represent  the 
characteristics  of  the  hybrid  processor.  For  example.  Table  8  shows  the  results  obtained  when 
using  a  sinusoidal  (i.e.,  multi-grey  level)  input  image  and  a  sinusoidal  error  function  for  the 
operatoi.  Again  convergence  of  the  aigoriihm  is  obtained.  In  particular,  we  see  achievement  of 
better  than  8-bit  accuracy  in  the  final  solution  after  only  four  error  correcting  iterations. 
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3.  Nonlinear  Optical  Devices  and  Digital  Optical  Computing 

^.1.  Si/PLZT  Opto-Electronic  Devices  for  POP 

Using  low  pressure  chemical  vapor  deposition  (LPCVD)  techniques  to  deposit 
polycTystalline  Si  onto  the  electro-optic  PLZT  substrate  and  then  recrystallize  the  polysilicon 
using  an  Ar'*'  laser  to  produce  single  crystal  silicon  grains,  we  have  fabricated  a  12  x  12 
electrically-addressed  spatial  light  modulator  (SLM)  array  and  a  16  x  16  optically-addressed 
SLM  array  [17],  The  schematic  diagrams  of  the  electronic  circuits  for  the  two  devices  are 
shown  in  Fig.  10. 

Each  pixel  of  the  ESLM  array  consists  of  a  driving  transistor,  a  reset  line  and  the  E-0 
modulator.  An  electrical  input  drives  the  NMOS  transistor,  the  output  of  which  in  turn  controls 
the  modulator.  In  experiments  performed  in  the  last  year,  we  have  achieved  excellent  transistor 
performance  characteristics  which  are  shown  in  Fig.  11  [18].  Figure  11a  shows  the  I-V 
characteristics  of  the  NMOS  transistor  and  Fig.  11b  shows  the  time  response  characteristics  of 
the  ESLM.  We  have  achieved  breakdown  voltages  of  48V  which  is  at  least  a  factor  of  2  greater 
than  that  of  previously  reported  thin  film  transistors  [19].  The  transconductance  was  200fJ.s 
which  then  yielded  a  calculated  value  of  550  cm^A^s  for  the  mobility.  We  are  at  present  in  the 
process  of  measuring  the  storage  time,  the  leakage  current  and  the  optical  modulation  depth  of 
these  modulator  arrays. 

The  performance  characteristics  of  the  OSLM  array  are  shown  in  Fig.  12.  In  this  device  we 
have  demonstrated  the  fabrication  of  a  CMOS  device  on  an  electro-optic  substrate.  Breakdown 
voltages  of  28v  were  achieved  for  the  NMOSFET.  The  transconductance  was  100^  sA^  and  the 
mobility  was  442  cm^A^s  [19].  The  major  improvement  of  the  performance  characteristics 
resulted  from  the  optimization  of  the  process  steps  and  the  improved  cry'Stal  grain  growth 
obtained  using  laser  recrystallization. 

3.2.  Digital  Optical  Processing 

To  study  optical  parallel  computing  architectures,  we  investigated  different  interconnection 
topologies  and  have  generated  holographic  optical  elements  (HOEs)  that  implement  (a)  fully 
interconnected  and  (b)  hypercube  interconnected  topologies  [20].  These  HOEs  were  space- 
variant  holograms  designed  and  generated  by  the  computerized  system  of  Fig.  1,  and  were 
recorded  in  dichromatic  gelatin  to  provide  diffraction  efficiency  of  better  than  80%. 
Reconstructions  from  one  sub-hologram  (PE)  for  fully  connected  and  one  for  a  hypercube 
interconnected  architecture  are  shown  in  Figs.  13  and  14,  respectively.  Recently,  we  are 
investigating  the  replacement  of  optical  holograms  by  computer  generated  holograms  for 
programmable  optical  interconnects  and  volume  holograms  in  photorefractive  crystals  for 
reprogrammable  optical  interconnects. 
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This  table  shows  the  mean  square  error  of  estimate  image  for  three  different  spike  removal  algorithms  and  seven  different  noise 
levels.  It  is  clearly  shown  that  the  performance  of  the  space-variant  filter  is  much  better  than  the  others,  especially,  for  the  very  noisy 
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Table  7. 

Comparison  of  results  for  the  BER  algorithm  obtained  from  a  theoretical  description  of  the 
algorithm  and  from  a  computer  simulation  of  a  hybrid  (electronic/optical)  implementation  of  the 
algorithm.  The  initial  input  was  a  constant  intensity  two-dimensional  image;  relative  operator  error 
(25%)  was  constant  over  the  input  image  domain.  Experimentally  realizable  results  will  be  limited 
by  the  accuracy  of  the  electronic  processor. 
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Table  8. 

Results  from  a  computer  simulation  of  a  hybrid  (electronic/optical)  implementation  of  the  BER 
algorithm  using  a  sinusoidal  input  image  function  and  a  sinusoidal  relative  operator  error  function 
(max  25%). 
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Fig.  1.  Schematic  of  the  hybrid  system  for  producing  computer  generated  space-variant 
holographic  filters.  A  matrix  of  256  by  256  holograms  of  different  PSF  is  recorded 
with  this  system.  T^e  following  focal  length  (f)  in  mm  and  F-number  (F#)  lenses  are 
employed  =  100,  fi  =  20.  /a  =  360,  /4  =  80  of  F#  =  1 .6,  /s  =  80  of  F#  =  1 .6,  fs  =  80 
of  F#  ^  1.6, /7  =  20, /s  =  120, /9  =  120. 


Fig.  2: 


on?he^Rmrr''r""H  ^ CTop  view).  TTie  input  pattern  is  displaved 
Dlaced  m  and  imaged  on  the  space-variant  holographic  filter.  A  camera  is 

placed  in  the  reconstruction  plane  of  the  holograms. 
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Fig.  5.  fn  r -6  coordinate  transform  of  a  rotation  variant  pattcm.  a)  Original  input  image  of 
an  off-center  circle,  b)  Optical  coordinate  transform  inr  -6  of  the  input  pattern,  c) 
Optical  coordinate  transform  of  the  input  pancm  rotated  by  an  angle  of  45°  around  the 
origin  of  the  coordinate  system,  d)  Optical  coordinate  transform  of  the  input  pattern 
rotated  by  an  angle  of  90°  around  the  origin  of  the  coordinate  system. 
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fnr-d  coordinate  transform  of  another  scale  variant  pattern,  a)  Input  image  of  an 
airplane,  with  original  scale  1.  b)  Optical  coordinate  transform  of  the  image  in  (a),  c) 
Input  image  of  (a)  scaled  by  a  factor  of  a  =  0.7.  d)  Optical  coordinate  transform  of  the 
image  in  (c).  The  transformed  pattern  of  (b)  is  shifted  upward  by  an  amount  of  fn  a. 


Fig.  6. 


(c)  (d) 


Fig.  7.  (a)  Original  noisy  image,  (b)  output  of  a  low-pass  filter,  (c)  output  of  a  median  filter, 

(d)  output  of  the  proposed  space-variant  filter. 
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Fig.  8a.  The  transmitted  output  signal  just  after  initial  illumination  of  the  stationary  input  scene. 

After  a  few  seconds,  the  output  vanishes  completely.  Notice  the  model  car  in  the  lower 
right. 


Fig.  8b.  The  transmitted  output  just  after  the  model  car  was  moved  to  the  left.  The  car’s  initial 
and  current  positions  are  revealed,  along  with  a  fading  trail  of  the  most  recent  posi¬ 
tions.  The  stationary  components  of  the  scene  remain  dark. 


(b) 


Fig.  9.  (a)  Tensor  storage.  The  Imth  component  of  W(l,m,y,y)  tensor  is  displayed  on  the 

electronically  addressed  spatial  light  modulator  (ESLM).  Lense  LI  images  W(l,m,.x.y) 
onto  the  random  phase  function  exp  lj0(,r,y)].  This  random  phase  function  is  provided 
by  a  high  resolution  E-Beam  fabricated  random  phase  code  transparency  (RPC).  TTie 
resulting  product  is  Fourier  transformed  via  lens  L2  onto  the  photorefractive  cr>’stal 
while  a  reference  plane  wave,  deflected  by  an  angle  la  and  ma,  arrives  from  the 
acousto-optic  device.  The  superposition  of  the  first  beam  and  the  shifted  reference 
beam  creates  the  appropriate  interference  hologram.  The  space  bandwidth  product  of 
the  photorefractive  crystal  is  high  (10^)  and  will  allow  the  storage  of  a  large  W(l.m,x,y) 
tensor,  (b)  Random  Phase  Coded  Matrix  Tensor  Multiplication:  When  G(x,y)  is 
presented  during  readout,  the  Im  correlations  of  W(l,mpc,y)  with  G(x,y)  sharpened  by 
the  random  phase  function  are  transmitted  along  the  various  (la,ma)  directions.  This 
beam  is  then  Fourier  transformed  by  lense  L3  and  captured  by  the  CCD  camera. 

After  the  initial  inputs  to  W(l,mpc,y)  are  made  the  two  phases  of  operation  will 
alternate.  An  image  will  be  input  and  the  correlations  determined.  Then  the 
appropriate  Imth  component  of  W  will  be  updated  by  an  amount  determined  by  the 
algorithm.  The  alternation  of  these  two  phases  provides  the  system  with  the  ability  to 
adjust  the  weighting  factors  and  to  learn  new  categorizations  as  they  appear. 


10b.  Schematic  of  the  fabricated  OSLM  circuit  in  each  unit  cell. 


Fig.  11a.  I-V  characteristics  of  a  typical  NMOS  pass  transistor  used  in  the  ESLM  circuit. 


Fig.  lib. 


Response  time  characteristics  of  ESLM.  The  upper  curve  is  plotted  on  a  longer  time 
scale,  the  lower  curve  shows  the  rise  and  fall  times  in  detail  on  a  shorter  lime  scale. 


Fig.  12  Typical  I-V  characteristics  of  Si/PLZT  NMOS  transistor  utilized  in  OSLM  array. 


Fig.  13.  Optical  reconstruction  from  one  of  the  32  x  32  arrays  of  holograms  for  full 
interconnection  among  the  32  x  32  processing  elements. 


